■HHHHBMMmm 

NPS  ARCHIVE 
1965 
LARSEN,  J. 


WW 


I 


"  MI|W" 


POLYCRYSTALUNE  ALUMINUM-COPPER 


JOHN  T.  LARSEN 


AN  X-RAY  METHOD  FOR  FOLLOWING  THE 

PRECIPITATION  PROCESS 
IN  POLYCRYSTALLINE  ALUMINUM-COPPER 


John  T.  Larsen 


AN  X-RAY  METHOD  FOR  FOLLOWING  THE 

PRECIPITATION  PROCESS 
IN  POLYCRYSTALLINE  ALUMINUM-COPPER 


by 

John  T.  Larsen 
Lieutenant,  United  States  Navy- 


Submitted  in  partial  fulfillment  of 
the  requirements  for  the  degree  of 

MASTER  OF  SCIENCE 

IN 
MATERIAL  SCIENCE 

Tihited  States  Naval  Postgraduate  School 
Monterey,  California 

19  6  5 


Me.  5 


Library 

U.  S.  Narval  Post  aradu ate  Scnooi 

Monterey.  California 

AN  X-RAY  METHOD  FOR  FOLLOWING  THE 

PRECIPITATION  PROCESS 
IN  POLYCRYSTALLINE  ALUMINUM-COPPER 

by 

John  T.  Larsen 

This  work  is  accepted  as  fulfilling 

the  thesis  requirements  for  the  degree  o] 

MASTER  OF  SCIENCE 

IN 

MATERIAL  SCIENCE 

from  the 

United  States  Naval  Postgraduate  School 


ABSTRACT 
An  experimental  study  has  been  made  of  the  feasibility  of 
following  the  precipitation  process  in  a  polycrystalline  aluminum- 
2%%   copper  alloy.   The  intermediate  structures  formed  during  ageing 
have  previously  been  determined  by  single  crystal  methods  to  pass 
through  the  fol lowing  sequence: 

The  sequence  of  structures  was  followed  by  analysis  of  the  diffrac- 
tometer  line  profile  of  the  ©c^OO)   peaks.   The  experimental 
profiles  were  corrected  for  instrumental  broadening  by  the  Stokes 
method.   One-dimensional  Patterson  projections  of  the  electron 
density  distribution  in  the  alloys  were  synthesized  from  the 
corrected  profiles. 

The  complete  sequence  of  intermediate  precipitates  could  not 
be  realized  in  this  investigation.   However  it  was  possible  by  this 
method  to  distingish  between  the  G.P.  \\~\  and  Q   structiires  by 

observing  the  height  and  shift  of  the  Patterson  function  origin  peaks. 
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1.     Introdxiction. 

The  phenomenon  of  precipitation  hardening  holds  great  commer- 
cial interest.   It  is  a  method  of  enhancing  the  mechanical  properties 
of  a  metal  without  resorting  to  other  processes  such  as  work  hardening, 
dispersion  hardening  or  a  martensitic  type  transformation.   It  elim- 
inates the  concern  regarding  the  hardenability  of  a  metal.   Metal 
parts  may  also  be  formed  or  fabricated  in  a  soft  condition  and  then 
hardened  by  a  relatively  simple  ageing  treatment  while  retaining 
excellent  dimensional  stability.   One  example  of  the  importance  of 
precipitation  hardening  is  its  application  to  aluminum  alloys  which 
has  resulted  in  a  great  number  of  usable  high  strength  materials. 

The  mechanism  of  precipitation  hardening  is  not  at  all  simple. 
The  most  basic  concept  involves  the  formation  of  a  stable  phase 
dispersed  throughout  a  matrix.   This  does  not  include  dispersion 
hardening,  which  is  due  to  the  presence  of  a  second  phase  dispersed 
in  the  matrix.   However,  it  is  not  a  precipitation  process.   Even 
in  the  basic  concept  of  precipatation  there  are  constraints  imposed 
on  the  nucleation  and  growth  of  such  a  phase  because  it  is  being 
formed  in  a  solid  solution.   In  general,  the  process  is  much  more 
complicated  than  this.   It  varies  with  each  alloy  system  but  usually 
involves  the  precipitation  of  a  metastable  phase,  which  later  gives 
way  to  a  stable  one.   There  may  be  as  many  as  three  intermediate 
stages  before  the  stable  precipitate  is  formed.   The  intermediates 
may  or  may  not  be  true  phases,  depending  on  the  alloy  system. 

The  intermediate  stages  hold  great  interest.   In  many  alloy 
systems  these  intermediates  have  the  greatest  influence  on  mechanical 
properties,  and  the  formation  of  a  stable  precipitate  results  in  a 


decrease  in  hardness  and  strength. 

The  Al-Cu  alloy  represents  an  ideal  system  in  which  to  study 
the  age-hardening  process.   It  is  representative  of  many  systems  for 
it  displays  both  substitutional  and  displacement  disorder  during 
ageing.   The  copper  atoms  have  a  greater  scattering  factor  and  a 
smaller  radius  than  the  aluminum  atoms  (tch/^^  =  29/13  and  2.56  A 
compared  to  2.86A  ).   The  difference  in  atomic  size  results  in  the 
displacement  disorder  while  the  greater  scattering  factor  of  copper 
allows  its  effects  to  be  detected  by  x-ray  analysis. 

The  age-hardening  system  in  Al-Cu  has  been  studied  more  than 
any  other.   As  a  result  its  sequence  of  structures  is  known  in  consid- 
erable detail. 

The  mechanism  of  the  decomposition  of  the  quenced  supersatur- 
ated solid  solution  shows  four  different  stages  whose  relative  roles 
depend  strongly  on  the  annealing  temperature.   In  general  the  sequence 
may  be  represented  as  follows: 

g.r  m  —*  c.p.\£i  -*■  e'  -*  © 

Two  stages,  9  and$,  correspond  to  true  precipitates  of  composition 
CuAlfc ,  where  Q   is  the  equilibrium  phase.   The  G.P.  refers  to  a 
"Guinier-Preston  zone",  which  may  be  defined  as  a  region  rich  in 
solute  atoms  which  is  formed  on  certain  planes  in  the  solvent  matrix 
and  retains  what  is  essentially  the  structure  of  the  matrix.   G.P.^^ 
is  a  copper  rich  region  of  plate-like  shape  formed  on  {VOO}  planes 
of  the  aluminum  matrix  whose  maximum  size  is  a  few  atomic  planes 
thick  and  up  to  about  100  A  in  diameter.   G«P.  £."£]  displays  features 


lj.  If.  Silcock,  T.  J.  Heal  and  H.  K.  Hardy,  J.  Inst.  Metals  82, 
240-242  (1953). 


of  a  zone  and  a  true  phase.   As  these  platelets  reach  a  thickness 
of  about  30  A  they  become  unstable  and  dissolve.   They  may  be  considered 
as  having  a  lattice  similar  to  that  of  the  aluminum  matrix.   To  a  first 
approximation  it  is  merely  a  superlattice  of  aluminum,  tetragonal 
with  c/a  =  2  where  a  is  the  parameter  of  aluminum. 

The  above  structures  do  not  necessarily  evolve  from  one  another 
in  sequence,  but  two  (or  more)  may  occur  simultaneously  in  the  same 
specimen.   The  conditions  have  been  determined  in  which  the  three 

metastable  phases  and  the  equilibrium  phase  appear  as  a  function  of 

o 
time  and  temperature  for  different  concentrations  of  Cu  in  Al . 

A  variety  of  x-ray  techniques  have  been  used  in  the  numerous 
precipitation  investigations.   Single  crystals  have  generally  been 
employed  using  both  white  and  monochromatic  radiation.   Some  work 
has  been  done  using  normal  powder  methods.   The  small  angle  scatter 
method  has  been  applied  with  great  success  by  Guinier.^  The  tech- 
niques using  single  crystals  and  monochromatic  radiation  have 
advantages  of  ease  and  certainty  of  interpretation  in  terms  of  a 
reciprocal  lattice  plot,  but  require  sophisticated  equipment  and 
precise  experimental  technique.   Some  of  the  results  obtained  by 
the  above  methods  have  recently  been  confirmed  by  electron  micro- 
scopic photographs.  **  The  use  of  an  x-ray  diffractometer  in  the 
analysis  of  polycrystalline  specimens  is  attractive  because  of  the 
relatively  simple  procedures  involved  and  the  possible  extension  of 


2 Ibid. ,  pp.  243-45. 
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the  technique  to  precipitation  studies  in  commercial  alloys.   Hence, 
it  is  the  purpose  of  this  investigation  to  follow  the  ageing  process 
in  a  polycrystalline  Al-4$Cu  alloy  by  means  of  a  diffraction  line 
shape  analysis  and  Patterson  projection  of  the  electron  density 
distribution  within  the  alloy. 

The  data  necessary  to  form  a  Patterson  function  is  available 
from  a  diffraction  line  profile,  as  is  shown  below.   However,  a  true 
line  profile  must  first  be  obtained  by  correcting  the  observed  profile 
for  instrumental  broadening,  due  mainly  to  the  slight  spread  of  wave- 
lengths in  the  incident  radiation.   A  method  has  been  derived  by 
Stokes  using  the  convolution  theorem  of  mathematics  to  correct  the 
experimental  line  shape  without  any  prior  assumptions  concerning  the 
corrected  line  shape.    The  method  consists  essentially  of  obtaining 
the  Fourier  coefficients  for  the  broadened  and  reference  curves, 
operating  on  these  coefficients  (termed  unfolding),  and  then  synthe- 
sizing the  corrected  curve  from  the  unfolded  coefficients.   This 
technique  eliminated  broadening  due  to  Koc  doubling  as  well  as  other 
forms  of  broadening. 

It  has  been  shown  that  the  true  line  profile  is  proportional  to 
the  Fourier  transform  of  the  square  of  the  absolute  magnitude  of  the 
structure  factor,  ^F^il  •    The  structure  factor  r*V\\cft»  in  general  a 
complex  quantity,  represents  something  real  about  a  crystalline 
material  for  it  is  the  Fourier  transform  of  the  electron  density  in 


5 
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the  material.^  Experimentally  only  |f-h\-^l   is  determined  which  yields 
no  information  concerning  its  phase.   If  the  phase  relations  could  be 
determined  experimentally  there  would  be  no  problem  in  finding  the 
arrangement  of  atoms  in  any  crystal,  no  matter  how  complicated.   Of 
the  methods  proposed  to  circumvent  the  phase  problem,  the  Patterson 
function  has  emerged  as  the  most  powerful  and  adaptable.   This  function 
is  a  Fourier  series  whose  coefficients  are  IFhfc.ft.l  .    If  one  considers 
the  Fourier  coefficients  \F yNVJj>\   as  derived  by  forming  products  of 
each  Fourier  coefficient  FvvVJ.  with  its  complex  conjugate,  the  electron 
density  function  can  be  derived  directly  from  the  Patterson  function. 
The  Patterson  function  is  essentially  a  display  of  all  the  vector 
lengths  between  atoms  of  a  unit  cell.   All  the  vectors  of  one  length 
will  then  appear  as  a  peak  in  this  display.   Since  the  reverse  vectors 
are  also  displayed,  it  is  necessary  to  sort  out  the  duplicated  peaks 
in  order  to  arrive  at  the  electron  density  function. 

2.    Alloy  preparation. 

The  Al-4$Cu  (by  weight)  alloy  was  prepared  using  99.99$  ingot 
Al  and  99.94%  Cu  shot.   The  aluminum  was  first  melted  in  an  induction 


heated  graphite  crucible,  and  then  the  copper  shot  was  added.   The 
molten  metal  was  held  between  700°  and  720  C  for  45  minutes  to  ensure 
complete  solution  of  the  copper.   The  alloy  was  cast  into  a  flat 
rectangular  zirconia  boat,  preheated  to  700  C,  and  allowed  to  solidify 
at  room  temperature.   The  ingot  thus  formed  weighed  125  grams, 


7 
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measuring  approximately  4"  x  2^"  x  ■§". 

Microscopic  examination  of  the  ingot  revealed  a  large  grained 
structure  containing  intergranular  theta  phase  and  no  evidence  of 
undissolved  copper. 

The  ingot  was  homogenized  at  542  C  for  24  hours  and  furnace 
cooled.   Grain  refinement  was  accomplished  by  cold  rolling  and 
recrystallizing  at  600  C  for  one  hour.   The  amount  of  reduction 
per  pass  through  the  rolls  was  limited  to  about  20$  due  to  work 
hardening,  so  the  refinement  treatment  had  to  be  repeated  three 
times,  resulting  in  an  ingot  .17"  thick,  and  containing  equiaxed 
small  grains  with  finely  divided  theta  phase  dispersed  throughout. 

Eight  hardness  specimens  were  cut  from  the  ingot,  each  being 
about  1.75"  x  .5"  x  .17"  in  size.   A  hole  was  drilled  in  one  corner 
of  each  specimen  and  top  and  bottom  surfaces  were  smoothed  on  no.  0 
and  no.  3/0  metal lographic  paper. 

3.     Solution  treatment  and  age  hardening. 

The  eight  specimens  were  strung  on  a  wire  with  a  £"  separation 

to  ensure  rapid  heat  transfer  both  on  heating  and  quenching.   They 

e 
were  solution  treated  in  a  salt  bath  for  six  days  at  536  C,  which 

is  about  36  above  the  solid  solubility  of  4%   copper  in  aluminum. 

Following  the  solution  treatment  the  specimens  were  quenched  in 

room  temperature  water,  and  the  as  quenched  hardness  was  taken 

immediately  on  one  specimen.   The  surface  irregularities  caused  by 

the  salt  bath  on  the  remaining  specimens  were  removed  by  lightly 


M.  Hansen,  Constitution  of  Binary  Alloys,  p.  84, 


polishing  on  no.  0  and  no.  3/0  paper.   They  were  then  suspended  in 
an  oil  bath  at  165  C  for  ageing.   The  as  quenched  hardness  specimen 
was  aged  at  room  temperature. 

Itficroscopic  examination  of  the  as  quenched  hardness  specimen 
following  solution  treatment  revealed  a  small  amount  of  undissolved 
theta  phase  in  the  structure.   A  lattice  parameter  measurement  on 
the  oC(422)  diffraction  peak  of  this  specimen  indicated  that  about 
2#fc   of  the  copper  had  gone  into  solution.   The  small  amount  of  undis- 
solved theta  would  have  no  effect  other  than  to  reduce  the  peak 
hardness  obtained  during  ageing. 

4.    Diffraction  procedure. 

As  each  specimen  reached  its  respective  ageing  time,  it  was 
quenched  in  room  temperature  water  and  hardness  readings  were  taken 
immediately.   No  less  than  six  readings  were  taken  on  both  sides  to 
ensure  a  good  average.   The  surfaces  were  then  given  a  rough  electro- 
polish  to  remove  the  oxide  layer  and  to  ensure  that  the  surface  was 
as  nearly  a  true  representation  of  the  bulk  of  the  specimen  as  possible 
Chemical  macro-etching  was  not  used  due  to  its  preferential  removal  of 
copper  from  the  alloy. 

The  o6(200)  peak  of  each  specimen  was  recorded  using  a  Norelco 
x-ray  dif fractometer  with  Ni  filtered  CuKoL  radiation.   The  diffracted 
radiation  was  counted  by  a  xenon  proportional  tube  in  conjunction  with 
a  pulse  height  analyzer  set  to  pass  the  central  90%  of  the  pulse 
distribution. 

The  reference  peak  for  the  Stokes'  unfolding  was  obtained  from 
a  recrystallized  specimen  of  the  pure  ingot  aluminum.   It  was  possible 


to  get  a  well  resolved  kot  doublet  for  the  (200)  peak  from  this  specimen. 

5.  Treatment  of  data. 

The  intensity  of  each  diffraction  peak  was  recorded  as  a  function 
of  '£©.   Due  to  the  sharpness  of  the  reference  peak  it  was  necessary  to 
use  an  interval  of  .0125  26  between  points  to  describe  the  peak  by  a 
smooth  curve.   It  was  possible  to  use  a  point  interval  of  .025  2.9  to 
describe  the  broadened  peaks.   A  total  interval  of  1.8   was  used  on 
all  broadened  peaks  except  the  as  quenched  and  aged  one  day  specimens, 
where  it  was  possible  to  use  slightly  shorter  intervals  due  to  their 
narrowness.   The  maximum  value  of  the  diffraction  peak  was  centered 
in  the  total  interval  in  all  cases. 

The  tinfolding  and  Patterson  function  synthesis  calculations  were 
performed  by  the  U.  S.  Naval  Postgraduate  School's  CDC  1604  electronic 
computer.   The  diffraction  intensity  data  was  punched  onto  computer 
cards  for  input  to  a  Fortran  program  developed  by  Prof.  J.  R.  Clark, 
Dept.  of  Material  Science  and  Chemistry,  to  execute  the  above  calcu- 
lations. 

6.  Experimental  results. 

ageing  time  curve  for  the  specimens  aged  at  165  C 
is  shown  in  Fig.  1,  page  9.   The  specimen  aged  at  room  temperature  for 
65  days  attained  a  Rockwell  "F"  hardness  of  77.4. 

The  lattice  parameter  of  the  as  quenched  specimen  was  4. 044  A, 

which  corresponds  to  2£  weight  per  cent  copper  in  solution.   The 

e 
specimen  aged  at  165  C  for  one  and  two  days  retained  this  same  lattice 

parameter.   The  lattice  parameter  of  the  specimens  aged  for  five  days 
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e 
and  longer  was  4. 049 A  ,  which  is  the  same  as  that  of  pure  aluminum. 

The  as  quenched  specimen  retained  the  same  lattice  parameter  on  ageing 
at  room  temperature  for  65  days. 

The  Patterson  functions  for  the  artificially  and  naturally  aged 
specimens  are  shown  in  Appendix  I  on  both  a  small  and  large  scale. 
The  ordinate  is  P(w),  where  w  is  the  interatomic  vector  length.   The 
vector  between  any  two  atoms  in  the  crystal  always  has  its  tail  at 
the  origin  in  the  Patterson  function.   The  abcissa  then  represents 
the  magnitude  of  w  as  it  ranges  from  zero  to  "a",  the  lattice  param- 
eter of  aluminum,  in  increments  of  .Ola.   The  Patterson  functions 
displayed  in  Appendix  I  are  the  one-dimensional  projections  of  the 
electron  density  in  the  <\00'>p,.  direction. 

7.     Discussion  of  results. 

The  lattice  parameter  and  hardness  measurements  were  used  to 
confirm  the  structures  present  at  each  ageing  time.   The  lattice 
parameters  were  determined  with  only  sufficient  precision  to  detect 
any  changes  as  a  function  of  ageing  time.   The  structures  contained 
in  each  specimen  were  determined  by  comparing  the  experimental 
hardness  curve  with  the  hardness  curves  of  Silcock,  et  al.     It 
was  not  possible  to  obtain  the  complete  sequence  of  intermediate 
precipitates  on  ageing  at  165  C  since  only  2$$>   copper  went  into 
solution.   Initial  ageing  at  165  C  produces  G.P.  \£\      and  ©  simul- 
taneously, while  at  peak  hardness  (five  days)  0  predominates.   The 
specimen  fully  aged  at  room  temperature  contains  only  G.P.  \X\      since 


Silcock,  et  al. ,  op.  cit.  ,  p.  243. 
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the  precipitation  process  does  not  progress  beyond  this  point  in 
naturally  aged  specimens. 

Lattice  parameter  measurements  confirmed  that  the  above  struc- 
tures were  present.   A  constant  value  of  parameter  is  characteristic 

12 
of  the  G.P.  zones.     It  was  observed  that  the  parameter  remained 

constant  during  natural  ageing.   The  parameter  also  remained  constant 
through  two  days  ageing  at  165  C,  indicating  that  G.  P.V^-1  was  present 
during  this  time.   At  five  days  ageing  the  parameter  had  relaxed  to 
the  pure  aluminum  value.   Hence  9  became  the  predominant  intermediate 
at  this  time. 

When  comparing  the  Patterson  origin  peaks,  Appendix  I,  Figs.  4 
through  11,  it  is  observed  that  the  as  quenched  and  naturally  aged 
specimens  have  peaks  exactly  at  the  origin  while  those  of  the  arti- 
ficially aged  specimens  are  shifted  slightly.   Also  the  origin  peak 
for  the  naturally  aged  specimen  is  roughly  twice  the  magnitude  of 
all  the  others.   It  has  been  shown  that  the  Patterson  peak  height 
is  not  a  true  measure  of  the  product  of  the  number  of  electrons, 

~t£\^-1  »  contained  in  the  two  atoms  separated  by  the  vector  that  corres- 

13 
ponds  to  the  peak.     The  peak  height  is  directly  proportional  to 

Z?^z  but  inversely  proportional  to  the  sum  of  the  cross-sections  of 

the  electron  density  peaks  for  a  one— dimensional  Patterson  peak. 

Therefore  a  high  Patterson  peak  is  not  necessarily  due  to  two  heavy 

atoms.   The  height  also  depends  on  the  product  of  the  electron 


11  Ibid- 

12 
"A.  Guinier,  Solid  State  Physics,  Vol.  9,  p.  341 

13Buerger,  op.  cit.  ,  p.  26. 
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densities  of  the  atoms  in  the  pair. 

The  copper  atom  has  a  greater  electron  density  than  aluminum. 
Copper  contains  29  electrons  compared  to  13  in  aluminum  and  is  11$ 
smaller  in  atomic  size.   Hence  the  presence  of  copper  in  an  atom 
pair  would  substantially  increase  the  Patterson  peak  height  for  that 
pair.   In  a  specimen  fully  aged  at  room  temperature  a  large  fraction 
of  the  copper  atoms  are  assembled  in  the  G.P.  zones,  which  are  a  few 
atomic  layers  thick  and  are  on  the  average  less  than  100A  apart.14 
When  looking  perpendicular  to  the  zones  in  the  (lOO)^  direction 
there  would  be  a  certain  number  of  vectors  of  length  "a"  between 
atom  pairs  containing  one  or  two  copper  atoms.   There  would  be  a 
smaller  number  of  like  vectors  in  the  as  quenched  specimen,  which 
is  assumed  to  be  a  random  substitutional  solid  solution  of  copper 
atoms  in  the  aluminum  matrix.   The  number  of  such  vectors  would  also 
be  less  in  the  artificially  aged  specimens,  which  initially  contained 
G.P.^J   and  0  and  later  onlyQ.   This  can  best  be  shown  by  consid- 
ering  the  accepted  structure  of  ©,  Fig.  2,  page  13.   The  0  struc- 
ture is  tetragonal  with  parameters  a  =  4.04A  and  c  =  5.8A,  and  the 
axes  are  parallel  to  the  OOO^p^  directions. 1^  Only  one  of  the 
three  possible  orientations  of  B  is  shown  in  Fig.  2.   When  all  three 
orientations  are  considered  together,  only  two  produce  copper  atom 
separations  of  4. 04  A  in  the  <\0O>^  direction.   Where  G.P.^ZQ  and 
B  occur  together  the  copper  atoms  are  shared  between  them  and  the 

e 

total  number  available  to  share  in  atom  pairs  4.04  A  apart  is  less 


14 

A.  Guinier,  Solid  State  Physics,  p.  341 

l  r 
°Silcock,  et  al. ,  op  cit.  ,  p.  242 
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than  in  the  naturally  aged  specimen.   This  is  also  true  when  only  Q 
is  present.   Hence  one  would  expect  G.P.[IJ  zones  to  produce  higher 
Patterson  origin  peaks  than  the  other  intermediate  structures. 

As  seen  in  Appendix  I  the  left  origin  Patterson  peaks  occur  at 
w  values  slightly  less  than  "a"  for  all  artificially  aged  specimens. 
This  shift  is  real  and  is  not  due  to  an  error  in  calculation  or 
plotting  of  the  Patterson  function.   The  peak  is  at  ,98a  for  one  days  " 
ageing  and  at  ,99a  for  longer  ageing  times.   This  is  best  explained 
by  considering  the  matrix  distortion  due  to  each  of  the  intermediate 
precipitates.   When  a  copper  atom  is  substituted  for  an  aluminum 
atom  the  matrix  lattice  contracts  locally  due  to  the  smaller  atomic 
size  of  copper.   All  three  intermediate  precipitates  are  coherent 
with  the  aluminum  matrix.   This  coherency  between  matrix  and  preci- 
pitate lattices  produces  distortions  in  the  matrix  removed  from  the 
precipitate  interface.   The  extent  of  the  matrix  distortion  depends 
on  the  size  and  shape  of  the  precipitate  particle.   Thus  when  G.P.^l"\ 
zones  are  formed  in  their  thin  disc  shape  the  lattice  distortions 
are  small  and  do  not  extend  far  from  the  interface.   The  average  cell 
parameter  is  then  that  of  the  matrix  and  the  Patterson  peak  remains 
at  \  w  \   =   a. 

When  G.P.LZ^  is  present  the  lattice  distortions  are  more  extensive 
than  those  due  to  G.P.Q\"^  ,  since  the  G.P.[<Q  zones  are  thicker  and  have 
larger  diameters.   The  distortions  due  to  0  are  also  extensive  because 
of  the  §  c  parameter  of  5.8A.   Thus  when  G.P.ljZ^  and  8  are  present 
together  after  ageing  one  day  the  lattice  is  distorted  by  their 
combined  effects  and  the  Patterson  peak  occurs  at  \w\  =  .98a.   As 
the  G.P.i.2^  zones  dissolve  and  1^  grows  with  longer  ageing  times  the 

14 


lattice  distortion  decreases  and  is  due  to  ©  alone  at  peak  hardness 
of  five  days. 

In  Appendix  I,  Figs.  12  through  19,  are  shown  the  Patterson 
functions  on  a  large  scale.   In  general,  the  number  of  resolved  peaks 
increases  as  the  ageing  time  increases.   This  is  due  to  the  increasing 
number  of  harmonics  used  in  the  synthesis  of  the  Patterson  function. 
The  number  of  harmonics  used  in  each  case  is  determined  by  the  point 
at  which  the  coefficients  of  the  broadened  curve  have  decreased  to  a 
small  value  and  those  of  the  reference  curve  still  have  an  appreciable 
value.   The  use  of  higher  harmonics  beyond  this  point  results  in  a 
greater  uncertainty  in  the  coefficients  of  the  corrected  curve.     The 
intermediate  peaks  in  the  above  Patterson  functions  are  all  quite  small 
compared  to  the  origin  peaks.   The  resolution  presented  by  the  inter- 

c 

mediate  peaks  is  considerably  greater  than  the  1,5  to  2A expected 

o 
using  x-radiation  of  1,542A  wave  length.   Thus  the  fine  detail  is 

false  and  is  only  a  result  of  the  harmonic  analysis.   No  information 

concerning  atomic  positions  can  be  gained  from  it.   It  is  observed 

that  the  fine  details  of  the  Patterson  functions  for  five  days  ageing 

and  longer  are  essentially  the  same  indicating  that  the  line  profiles 

are  very  nearly  the  same. 

8.     Conclusions  and  recommendations. 

The  scattered  intensity  distribution  in  reciprocal  space  for 

17 
the  aluminum-copper  alloy  is  shown  in  Fig.  3.     The  scattering 


16Stokes,  op_.  cit. ,  p.  389. 

17 

A.  Guinier,  Solid  State  Physics,  p.  339. 
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streaks  around  the  nodes  of  the  reciprocal  lattice  are  due  to  the 
destruction  of  the  strictly  periodic  electron  density  distribution 
in  the  crystal.   This  destruction  is  caused  by  the  G,  P.  zones. 


002 


000 


200 


Fig.  3.   Scattering  in  the  reciprocal  space  of  Al-Cu. 

The  streaks  pass  precisely  through  the  nodes  of  the  matrix  and  along 
the<VOO>  rows  of  the  reciprocal  lattice.   The  fact  that  the  G.P. 
zones  distort  the  matrix  lattice  causes  the  scattering  streaks  to 
appear  only  on  the  high-angle  side  of  the  reciprocal  lattice  nodes. 

It  was  anticipated  that  the  diffuse  streaks  about  the  nodes 
would  contribute  to  the  diffraction  process  in  a  distinctive  way.   It 
should  be  possible  to  differentiate  G.P.O"}   from  G.P. ^21   since  the 
length  of  the  diffuse  streaks  are  related  to  the  size  of  the  zone. 
This  distinction  could  not  be  realized  in  the  present  investigation 
as  G.P.  [2.J  zones  were  not  obtained  by  themselves  on  ageing.  However 


16 


it  was  possible  to  distinguish  G.P.[J3   from  0  by  observing  the 
height  and  shift  of  the  Patterson  origin  peak. 

The  amount  contributed  to  the  diffraction  profile  by  the  diffuse 
scattering  streaks  is  uncertain.   It  is  possible  that  the  broadening 
of  the  profile  due  to  instrumental  factors,  which  were  corrected  for, 
is  greater  than  that  resulting  from  the  diffuse  streaks.   However, 
since  the  instrumental  broadening  remains  constant  with  ageing  the 
profile  changes  due  to  the  diffuse  scattering  streaks  could  be  detec- 
table.  The  thermal  vibrations  of  the  atoms  about  their  average  position 
in  the  crystal  also  produce  anisotropic  diffuse  diffraction  effects 
about  the  nodes  of  the  reciprocal  lattice.   The  thermal  diffuse 
scattering  is  related  to  the  elastic  properties  of  the  crystal,  which 
certainly  change  during  ageing.   Hence,  the  changes  in  thermal 
scattering  should  also  contribute  to  the  line  profile  changes  during 
ageing.   The  amount  contributed  and  the  best  means  of  detecting  it  is 
a  difficult  problem  to  analyze  quantitively.   There  has  been  no 
attempt  in  this  investigation  to  do  so. 

There  are  several  ways  in  which  the  application  of  this  method 

of  analysis  could  be  improved.   The  use  of  monochromatic  x-rays  would 

eliminate  the  profile  broadening  due  to  the  spectral  spread  of  the 

incident  beam.   In  order  to  observe  the  effect  of  each  intermediate 

precipitate  on  the  line  profile  it  is  necessary  to  have  them  existing 

separately  in  the  alloy.   The  transition  in  the  aluminum-copper  system 

from  G.P.^Q  — *  G.P.\j£^  could  best  be  observed  by  ageing  a  4%   copper 

alloy  at  110  C,  and  from  G.P.^^-1*  0  by  ageing  the  same  alloy  at  190  C. 

The  complete  intermediate  sequence  could  be  achieved  with  definite 
separations  in  structures  by  ageing  a  4.5%  copper  alloy  at  130  C  but 
would  require  ageing  times  of  about  400  days. 

17 
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APPENDIX  I 
PATTERSON  PROJECTIONS 

The  Patterson  functions  for  the  naturally  and  artificially 
aged  specimens  are  shown  on  a  small  scale  in  Figs.  4  through  11, 
and  on  a  large  scale  in  Figs.  12  through  19.   The  ordinate  is  P(w), 
where  w  is  the  interatomic  vector.   The  abcissa  represents  the 
magnitude  of  w  as  it  ranges  from  0  to  "a". 
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APPENDIX  II 

FOURIER  COEFFICIENTS  FOR  REFERENCE,  BROADENED 
AND  CORRECTED  PROFILES 


The  column  headings  for  the  tabulated  Fourier  coefficients 
are  as  follows: 

N,  the  harmonic  number 

ALPHA,  coefficient  of  the  real  component 
BETA,  coefficient  of  the  imaginary  component 
INTENSITY  for  each  harmonic  is  computed  by 

where  L  is  the  number  of  equal  intervals 
between  the  data  points  used  to  describe 
the  curve. 
The  corrected  coefficients  are  normalized.   The  arrow  in  the 

"corrected"  column  indicates  the  number  of  harmonics  used  in  the 

synthesis  of  the  Patterson  function. 
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APPENDIX  III 
EXPERIMENTAL  DIFFRACTION  LINE  PROFILES 

The  following  instrumental  settings  were  used  in  obtaining  all 
of  the  diffraction  line  profiles.   The  abcissa  is  in  degrees  20. 

Cu  K<*  radiation,   50  KV,   40  ma,   Ni  filter 

•  9 
Scan  Speed  \  /minute,   Chart  Speed   30"/hour 

Time  Constant   2  seconds,    Scale   64  x  .6 

o 
Angular  Aperture   1  ,     Receiving  Slit   0.003" 

Xenon  Proportional  Counter,   1725  volts 

Base   7  volts,    Window   3  volts 
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THE  FOLLOWING  IS 
THAT  WILL  PERFORM  THE 
FUNCTION  SYNTHESIS. 


A  FORTRAN  60  COMPUTER  PROGRAM 
STOKES  UNFOLDING  AND  PATTERSON 
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OF  A  FOURIE 
SULTS  ARE  PR 
IMENSIONAL  S 
BULAR  FORM 
IMENSIONAL  S 
PERIMENTAL  V 
IZATION*  MUL 


COMPUT 
OEFFIC 
OEFFIC 
OEFFIC 
OURIER 
GRAVI 
OF  THE 
URVES 
BE  PRE 
ENTS  I 
N.  Y(U 
R  TRAN 
ESENTE 
ECTION 


E 

IENTS  A(N,1)  AND  B(N,1)  OF  A  CURVE 

IENTS  A(N»2)  AND  6(N,2)  OF  A  SECOND  CURVE 

IENTS  A(N)  AND  B(N)  OF  THE  UNFOLDED  CURVE  2/1 

TRANSFORM  (SYNTHESIS)  OF  THE  UNFOLDED  CURVE 
TY  OF  THE  UNFOLDED  CURVE 

MAXIMUM  IN  THE  UNFOLDED  CURVE 
CAN  BE  SUCCESIVELY  UNFOLDED 
SENTED  IN  BOTH  TABULAR  AND  GRAPHICAL  FORM 
F(H»K.»L)I**2»  THE  PROGRAM  COMPUTES  THE 
tV.W)  (I.E.  THE  SUMMATION  OVER  THE  COSINE 
SFORM  OF  I (F, (H,K,L) 1**2  >• 
D  AS  ONE  OR  TWO  DIMENSIONAL  SECTIONS 
S  CAN  BE  PRESENTED  IN  BOTH  GRAPHICAL  AND 


ECTIONS  ARE  PRESENT 
ALUES»  IF(H,K»L)I** 
TIPLICITY  AND  TEMPE 


ED  IN  TABULAR  FORM 
2  ARE  CORRECTED  FOR 
RATURE  FACTORS 


LORENTZ, 


FIRST  DATA  CARD  CARRIES  GENERAL  TITLE  -  COLUMNS  2-80 


SECOND  DATA 


CARD 
COL 


COL  4-5   NB  = 


1-3;  NA  ='  1  FOR  THE  FIRST  CURVE  OF  A  SET 
=  2  FOR  THE  SECOND  CURVE  OF  A  SET 

0   IF  FOLDING  IS  TO  BE  CALCULATED 
=  i  TO  COMPUTE  FOURIER  COEFFICENTS  ONLY 
A  VALUE  FOR  NB  MUST  BE  ENTERED  OR  STOP  15 
WILL  RESULT. 

DIMENSIONS  OF  AA,  BB,  AND  FIX  DEPEND  ON 
VALUE  OF  NA  SPECIFIED 
=  +1   SYNTHESIS  OF  UNFOLDED  DATA 
=  -1   SYNTHESIS  OF  NEW  DATA 


COL  6-7   NC 
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COL  8-9   ND  =  0    NO  GRAPH  OF  PATTERSON  FN. 

=  1   GRAPH  PATTERSON  FN.   SEE  BELOW. 
COL  10-74     TITLE  OF  UNBROADENED  CURVE 

THIRD  DATA  CARD  -  CONSTANTS  FOR  ANALYSIS  OF  UNBROADENED  CURVE 

COL   1  -   5   N  =  NO  OF  INTERVALS 

6-10   NCG  =  -1  FOURIER  SYNTHESIS  -  GRAPH 
=   0  NO  SYNTHESIS 
=  +1  FOURIER  SYNTHESIS 
11  -  15   MODE  =  1  IMPUT  DATA  PRINTED 

=  0  IMPUT  DATA  NOT  PRINTED 
16-20   J  -  FIRST  OUTPUT  HARMONIC 

(NORMALLY  -  0) 
21-25   K  -  HARMONIC  INCREMENT 
(NORMALLY  =  1) 
.  26-30   L-  FINAL  HARMONIC 

31  -  35   K2-  SUB  HARMONIC  -  IF  0»  PROGRAM 

SETS  K2  =  1  . 
36  -  40  '  LC  =  0   READS  DATA  CARD 
MEMORY  - 
=  1   USES  DATA  AVAILABLE  IN 
41  -  50   THO  =INITIAL  TWO  THETA- 
51  -  60   THF  =  FINAL  TWO  THETA 
61  -  73   XL  =  WAVELENGTH 

FOURTH  DATA  CARD  -  UNBROADENED  CURVE»N  EQUALLY  SPACED  POINTS  FORMAT 
(16F5.1)  ,  ON  THE  NECESSARY  NUMBER  OF  CARDS. 

REPEAT  SEQUENCE  WITH  A  DUPLICTE  GENERAL  TITLE  CARD  »  A  SIXTH 
CARD  SETTING  NA  =  2  AND  CARRYING  A  TITLE  FOR  THE  BROADENED 
CURVE  t  A  SEVENTH  CARD  SETTING  THE  CONSTANTS  FOR  ANALYSIS  OF  THE 
BROADENED  CURVE  »AND  AN  EIGHTH  CARD  WITH   THE  N  POINTS  OF  THE 
BROADENED  CURVE. 

FOR  A  GRAPH  OF  THE  PATTERSON  FN,  FOLLOW  THE  TWO  SETS  OF  DATA  CARDS 
(FOR  THE  REFERENCE  AND  BROADENED  CURVES)  WITH   TWO  MORE  CARDS. 
THE  FIRST  ONE  IS  BLANK  AND  THE  SECOND  ONE  SETS  NA=NB=NC=0t   ND=1. 

•  FURTHER  CONVOLUTIONS  CAN  BE  CARRIED  OUT  WITH 
ADDITIONAL  SETS  OF  FOUR  DATA  CARDS  »  INCREASING  NA  BY  ONE  FOR 
EACH  FOLD  UP  TO  10. 

ADDITIONAL  CURVES  CAN  BE  UNFOLDED  BY  ADDING  ADDITIONAL  SETS  OF 
DATA  CARDS  SETTING  NA  =1  AND  2  IN  THE  SECOND  AND  SIXTH  CARDS  OF 
EACH  SET 

FOR  A  FOURIER  ANALYSIS  SET   NB  =  1  IN  COL  5  IN  THE 
SECOND  CARD  OF  A  SET  OF  DATA  CARDS 

\ 


c 

C  THE  PROGRAM  IS  TERMINATED  BY  TWO  BLANK  DATA  CARDS 

C 

C 

C  STATEMENT  NUMBERS 

C    1  2    3    4    5    6    7    8    9   10   1 1   12   13   14   1 5   16   17 

C   18  19   20   21   22   23   24   25   26   27   28   29   30   31   32   33   34 

C   35  36   37   38   39   40   41   42   43   44   45   46 

C  51   52   53                                           / 

C   60  61 

C   <  103  104  105  406  107  108  109 

C  110  111  112  113  114  115  116  117  118  119 

C  120  121  122  123  124  125  126  127  128  129 

C  130  131  132  133  134  135  136  137  138  139 

C  140  141  142  143  144  145  146  147  148  149 

C  150  151  152  153  154  155  156  157  158  159 

C  160  161  162  163  164  165  166  167 

C  501  502  503  504  505  506  507  508  509  510  511  512  513  514  515 

C  900  901  902      904      906  907  908  909 

C  47  48  49  70  170  71  72  73 

C 

C 

C 

1  FORMAT  (1H1////) 

2  FORMAT  (80H0 

1     i  ) 

24  FORMAT     (I5»    2F10.4,    F10.5     ♦    2F10.4    ,F10.5    ,2F10.4,    F10.5) 

39  FORMATU15,   2F15.4,  F15.4//) 

40  FORMAT( 1H013X1HN10X5HALPHA11X4HBETA6X  9HINTENSITY) 

900  FORMAT   (  13,  3  12,  8A8 ) 

901  FORMAT(8I5.  3F10.3)  • 

902  FORMAT  (16F5.1) 

904  FORMAT(3H0N=I4,   7H      J=I5»   7H      K=I5»   7H      L=I5» 
1   8H      K2=I5»7X6HNCG  =  I3///12H  DATA  POINTS//) 

906  FORMAT  ( 1H0  16X  9HREFERENCE  20X  11HBROADENED    21X  9HCORRECTED  //) 

907  FORMATdHO  4X  1HN  2X  5HALPHA  7X  4HBETA  3X  9HINTENSITY  4X  5HALPHA 
1  ,6X  4HBETA  2X  9HINTENSITY  4X  5HALPHA  6X  4HBETA  2X  9H INTENS I TY// ) 

908  FORMAT (IX) 

909  FORMAT( IX,  16F5.1) 


C 

c 


DIMENSION   Y(200),  AA(200,4),  BB(200,4),  FIX(200,4),  NAME(8), 

1  FR(200,4),   FI(200,4),   FIT(200,4), 

2  Z(l, 1,900) 


51 


DIMENSION 
DO  71  1=1 
71  ITITLE(I)  = 
ITITLE(l)  = 
ITITLE(2)  = 
ITITLE(7)  = 
ITITLE(8)  = 
ITITLE(9)  = 
ITITLE(IO)  = 


WW(900) , 

12 
8H 

8H  LARSEN 
8H   0278 
8HPATTERSO 
8HN   FUNCT 
8HION   FOR 
8H   AL-CU 


ITI*TLE(12) 


/ 


LABEL  =  4H  65D 


47 

I  15 

49 

i  48 
9 

:  31 
30 


33 
32 


12 


515 
513 


511 

514 

512 

13 


MODE.  Jt  K»  L»  K2,  LC ♦  THO,  THF,  XL 


TPI  =  6.28318530718 

R=l.  0/(57. 2957795 13 1*2.0) 

READ  2 

READ  900,  NA»  NB,  NO  ND,  NAME 

IF(NC)  122,  15,  47 

IF(NA)  15.  15,  9 

■IF  (ND)  48,48,49 

CALL  DRAW(  101  ,Z  ,WW  ,0  ,0  ,  LABEL,  I  T  ITLE  ,0  , 1  0.0  ,0  ,0  ,0  ,0  ,9  » 10.,  1  ,LAST  ) 

GO  TO  3 

STOP  48 

READ  901,  N,  NCG, 

IF(K2)  30,31,30 

K2  =  1 

FN  =  N 

TNR  =  2. /FN 

THF  =  THF*R 

THO  =  THO  *  R 

IF  LC  EQUALS  ZERO, 

IF  (LC)  32,33,32 

READ 

PRINT  1 

PRINT  2 

PRINT  900,  NAi  NB, 

IF  (MODE)  12,13,12 

PRINT  904,   N,  J,  K,  L,  K2 , 

PRINT  909,  (Y( I ) ,  I  =  1,N) 

IF  (NA  -  1)  515,  513,  514 

STOP  515 

AREA  =  0.0 

DO  511  I  =  1,N 

|AREA1  =(Y(  I  )*(SINF(THF)  -  S I NF ( THO ) ) ) / 

AREA  =  AREA  +  AREA1 

DO  512  I  =  1,N 

Y(  I  )  =  Y( I )  /  AREA 

XK2  =  K2 

XJ  =  J 

ARG  =  TPI/FN/XK2 

ARG1  =  ARG*XJ  a 


READ 
902, 


NEW  DATA,  OTHERWISE  OPERATE  ON  OLD 
( Y( I ) ,  1=  1,N  ) 


NC,  ND, 


NAME 
NCG 


(XL  *  FN) 
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CS  =  COSF(ARGl) 

TCS  =  2.*CS 

SS  =  SINF(ARGl) 

XK  =  K 

ARG2  =  ARG*XK 

CB  =  COSF(ARG2) 

S(B  =  SINF(ARG2) 

M1=J 

GO  TO  (501*502, 503.504, 505,506,507,508, 509, 510) »  NA 

501  NAN=NA 
GO  TO  14 

502  NCW=NA 
MAN=NA-1 
GO  TO  14 

503  N;CW  =  NA-1 
NAN=NA-2 
GO  TO  14 

504  NCW=NA-2 
NAN=NA-3 

GO  TO  14  N 

505  NCW=NA-3 
NAN=NA-4 
GO  TO  14 

506  NCW=NA-4 
NAN=NA-5 
GO  TO  14 

507  NCW=NA-5 
NAN=NA-6 
GO  TO  14 

508  NCW=NA-6 
NAN=NA-7 
GO  TO  14 

509  NCW=NA-7 
NAN=NA-8 
GO  TO  14 

510  NCW=NA-8 
NAN=NA-9 

BEGINNING  OF  LOOP  ON  Ml 
14  T3  =  0. 
T2  =  0. 
I  =   N 
.  18  Tl  =  TCS.  *  T2  -  T3  +  Y(  I  ) 
I  =  1-1 
IF  (I)  16,19,17 

16  STOP  16 

17  T3  =  T2 
T2  =  Tl 
GO  TO  18 

19  M=M1+1 


5^ 


i 


/ 


AA(M,NA)  =  (Tl*-  CS*T2)*TNR 

BB(M»NA)  =  (SS*  T2)*TNR 

FIX(M,NA)  =  (  AA(iM,NA)*AA(M,NA)  +  BB(M,NA)  *  BB ( M,NA ) ) /TNR 

TEST  FOR  END  OF  Ml  SERIES 

IF(Ml-L)  21.26,20 

20  STOP  20 

21  M1=M1+K 
TCS  =  CB*CS  -  SB*SS  ' 
SS  =  CB*SS  +  SB*CS 
CS  =  TCS 
TCS  =  TCS  +  CS 
GO  TO  14 

26  NO=NA-l 
11  IF(NO)  23,42,  27 
23  STOP  23 

42  IF(NB)  27,3,27 

27  PRINT  1 
PRINT  2 
L2=L+1 
J2=J+1 
IF(NB)  37,38,37 

37  PRINT  40 
PRINT  908 
DO  43   M  ■  J2,  L2 
M2=M-1 

43  PRINT  39»M2,AA(M,NA) »BB(M»NA) »FIX(M,NA) 
QO  TO  3 

38  PRINT  906 
PRINT  907 
IF(NA-2)  22,22,28 

28  AA(M,NAN)=AA(M,NO) 
BB(M,NAN)=BB(M,NO) 
FIX(MtNAN)  =  ( AA(M,NAN)*AA(M,NAN)+  BB < M,NAN ) *B8 ( M»NAN ) ) /TNR 

22  DO  25  M=J2,L2 
M2=M-1 
J2  =  J  +  1 
FR(M,NO)  =  (AA(M,NCW)* 

1  *TNR 

FI(M»NO)  =  (BB(M,NCW)* 
1  *TNR 

IF(M  -J2)46,  45,  46 

45  FA  =  FR(J2,NO) 

46  FR(M,NO)  =  FR(M,NO)/FA 
FI (M,NO)  =  FI (M,NO)/FA 

FIT(MfNO)  =(FR(M,NO)*FR(M,NO)  +  FI(M,NO)*  F I ( M, NO ) ) /TNR 
PRINT  24  ,M2,AA(M,NAN)  ,BB(  M,NAN  )' ,F  I  X  (  M,  NAN  )  ,  AA(M,NCW),  BB(M,NCW 
1)  ,  FIX(M,NCW)  ,FR  (M,NO)  ,  FI(M,NO)  ,FIT(M,NO) 
IF(XMODF(M  ,45) )29,35,29 
35  PRINT  1 


AA(M,NAN)+  BB(M,NCW)*6B(M,NAN) ) /FIX (M, NAN) 
AA(M,NAN)-  AA(M,NCW)*  BB ( M, NAN ) ) /F I  X ( M, NAN ) 


GO 


c 
c 
c 
c 
c 
c 
c 
c 
c 
c 
c 
c 
c 
c 
c 
c 
c 
c 
c 
c 
c 
c 
c 
c 
c 
c 
c 
c 
c 
c 
c 
c 
c 
c 
c 
c 
c 
c 
c 
c 


PRINT 
PRINT 
PRINT 
GO  TO 


2 

906 
907 
25 


29  IF(XMODF(M 
36  PRINT  908 

GO  TO  25 
25  CONTINUE 

IF(NCG)  122 


>5)  )25»36»25 


44,  122 


PROGRAM  FOUSYN 


1ST  DATA  CARD 

.  ;COL  1-80  TITLE 

2ND  DATA  CARD 

COL  1-4  DEBYE  TEMP.  B 

B=0.0  IF  TEMP, 
11  LAMBDA 


FACTOR  IS  NOT  COMPUTED 


COL 


5  - 

11 

15 

21 

26 


-  14  NUMBER  OF  PEAKS  USED 

-  20  NUMBER  OF  INTERVALS/CELL  IN  X 

-  25  NUMBER  OF  INTERVALS/CELL  IN  Y 

-  30  NUMBER  OF  INTERVALS/CELL  IN  Z 


DIRECTIONS 
DIRECTION 

DIRECTION 


ADJUST  DIMENSION  STATEMENT  FOR  FRACTIONS  UVW  OF  XYZ  USED 


3RD  DATA  CARD 

COL  1-3  MAX  VALUE 

4-6  MAX  VALUE 

MAX  VALUE  OD  W 


OF 
OF 


7-9 


10 
13 
16 
19 
21 
24 


12 
15 
18 
21 
23 
26 


MAX 
MAX 
MAX 
MAX 
MAX 
MAX 


VALUE 
VALUE 
VALUE 
VALUE 
VALUE 
VALUE 


U 
V 

OF 
OF 
OF 
OF 
OF 
OF 


FOR 
FOR 
FOR 


PRINTOUT 
PRINTOUT- 
PRINTOUT 


ACCORDING  TO  SECTIONS 


WANTED 


ADJUST  DIMENSION  STATEMENT  FOR  MAX  VALUES  OF  HKL  USED 


4TH  DATA  CARD 
COL  1  -  3  H 
4  -  6  K 
7  -  9  L 
10  -  14 
15-18 
19  -  20 


RELATIVE  INTENSITY  (WITH  DECIMAL  POINT) 

TWO  THETA 

MULTIPLICITY 


61 


c 
c 

c 

c 


122  CONTINUE 


12,  F5.0) ) 


51  FORMAT  (3(313,  F5.0,  F4.0, 

52  FORMAT(F4.2,F7.5,  14,  315  ) 

53  Format  (613,  3F4.2,  2F10.4)  . 

4  FORMAT  (1H0,2X,  1HN,  6X,  22H(    U   ,    V  '  ,    W 
1     //) 

5  FORMAT  (15,  6X,  F6.4  ,  IX,  F6.4,  IX,  F6.4,  6X,  F12.4) 

6  FORMAT(50H 

.  7  FORMAT  (IX) 

8  FORMAT  (1H1////  ) 

60  FORMAT(21HOCENTER  OF  GRAVITY  =  F6.2) 
61 


)9X,  9H  Y(U,V,W) 


FORMAT  (55H0MAX 
13  //) 


HARMONIC  USED  FOR  SYNTHESIS  OF  CORRECTED  CURVE  = 


167  DIMENSION 
1   S3(100), 


F(l, 1,500),  TTH(  1,1,500)  ,  MU  (  1  , 1,500),  C3(100.)» 
Fid,  1,500),    Cl(l),   C2(l),   SKI),   S2  (  1  ) 


166 
129 


128 


136 


IF(NC)  129,  166,  128 

STOP  166 

READ  6 

READ  52,B,XL 

READ  53,  Jit 

GO  TO  136 

NN1  =  0 

NN2  =  0 

NU  =   0 

NV  =   0 

NW  =    N 

Jl  =  0 

J2  =  0 

J3  =  100 

Jll  =  Jl  +1 

J22  =  J2  +1 

J33  =  J3  +1 

M2  =  M  -  1 

UMAX  =  0 

•I2MAX  =  0 

I3MAX  =  M2 

Ul  =  0.0 

VI  =  0.0 

Wl  =  FN 

TTHO  =  THO 

TTHF  =  THF 

I1M  =  UMAX  +  1 

I2M=I2MAX+1 


,NP,NU,NV,NW 
J2»  J3»  UMAX,  I2MAX,  U3MAX,  Ul ,  VI,  Wl,  TTHO,  TTHF 


6Z 


I-3M  =  I3MAX  +1 
TTHD  =  TTHO  -  TTHF 
DO  115  13=1, I3M 
IS3  =  (  13  -  1)*100 
DO'  115   12  =  It  I2M 
IS2  =  (  I2-l)*10  +  IS3 
DO  115  11  =  1,  I1M 

IS=IS2+I1  / 

F(IS)=0.0  / 

Fl ( IS)  =  0.0 
TTH(IS)=0.0 
MU  (  I S ) =  0 . 0 
115  Z( IS)=0.0 

IF!  (NO  12  7,   15,  126 

126  DO  137  M  =  1,  L2 

IF  (M-l )  159,159,160 

159  I4M    =0 

160  IFU4M)     161,161,158 

161  F( 1,  1,  M)  =  FR(M,  NO)  j 
Fl  (1,  1,  M)  =  FI  (M,  NO) 

IF(FR(M,  NO)  )  170,  170,  137 
170  I4M  =  M 

GO  TO  137 
158  F( 1,  1,  M)  =  FR(M,  NO) 

Fid,  1,  M)  =  FI  (M,  NO) 
137  CONTINUE 
C        FOR  TRY  =  -1.0,   Y(U,V,W)  IS  SYNTHESIZED  USING  THE  FIRST  M 
C        HARMONICS  THAT  KEEP  Y(U,V,W)  POSITIVE.' 

162  TRY  =  -1.0 
I3M  =  I4M  -1 
I5M  =  I4M  -2 
GO  TO  125 

C         FOR  TRY  =  1.0,    Y(U,V,W)  IS  SYNTHESIZED  USING  ALL   L  HARMONICS 

163  TRY  =  1.0 
I5M  =  L2  -  1 
I3M  =  L2 

GO  TO  125 

127  IF  (XMODF(NP,  3 )  ) 1 1 1 »  112»  111 

112  NP  =  NP/3 
GO  TO  113 

111   NP  =  (NP/3)  +  1 

113  DO  114  NR  =  It  NP 

READ  51,  NN1,  NM2,  NN3,  F(IS)t  TTH(IS),  MU  (  I S )  ♦  FKIS) 

IS=NN3*100+NN2*10+NN1+1 

TTH( IS)=TTH( IS)*R 

AAA=SINF(TTH( IS) ) 

AAA2=AAA**2 

TE=1.0/EXPF(B*AAA2/XL**2) 

ZZ=COSF(2.0*TTH( IS) ) 

G3 


114 
125 

131 

130 

132 


138 
139 


157 

140 

154 
155 

153 


POL=(  1.0+ZZ**2)/AAA2*COSFr(TTH(  IS)  ) 

AMU( IS)  =  MU( IS) 

Fl(IS)  =  F1(IS)/(TE  *  POL  *AMU(IS)) 

F(IS)  =  F(IS)/(TE  *  POL  *AMU(IS)) 

PRINT  8 

IF  (NO  131,  15,  130 

PRINT  6 

GO  TO  132 

PRINT  2 

PRINT  61,  I5M 

PRINT  4 

TPI=6. 28318530718 

MARK=0 

IF(NC)  157,  138,  139 

STOP  138 

JU  =  1 

UV  =  1 

GO  TO  140 

DO  104   JU  =  1,  Jll 

DO  104  JV  =  1,  J22  v 

DO  104  JW  =  1,  J33 


IF(NC)  153, 
STOP  154 
U  =  0.0 
V  =  0.0 
GO  TO  156 


154,  155 


AJU  =  JU  -  1 
ANU  =  NU 
U  =  AJU/ANU 
AJV  =  JV  -  1 
ANV  =  NV 
V  =  AJV/ANV 
156  AJW  =  JW  -  1 
ANW  =  J  3 
,W  =  AJW/ANW 
IF(NC)  141,  142,  143 

142  STOP  142 

143  II  =  1 
12  =  1 

GO  TO  144 
141  DO  103  11  =  1,  I1M 
DO  103  12  =  1, I2M 

144  DO  103  13=1, I3M 
NN1  =11-1 
NN2  =12-1 
NN3  =13-1 
FN1  =  NN1 

FN2  =  NN2 
FN3  =  NN3 


€>4- 


IF  (NO  145,  146,  147 

146  S/OP  146 

147  Sl( II)  =  1.0 
S2(  12)  =  1.0 
GO  TO  148 

145  SKID  =  SINF(U*FN1*TPI  ) 
S2(  12)  =  SINF(V*FN2*TPI ) 

148  S3(I3)  =  SINF(W*FN3*TPI )  / 
TS  =  SKID  *  S2(I2)  *  S3(I3)  *  FKI1,I2,I3)          / 
IF  (NO  149,  150,  151 

150  STOP  150 

151  C 1 ( 1 1 )  =  1.0 
C2  (  12)  =  1.0 
GO  TO  152 

149  CKID  =  C0SF(U*FN1*TPI  ) 
C2(I2)  =  C0SF(V*FN2*TPI ) 

152  C3(  I3)=COSF(W*FN3*TPI  ) 

T  =  CK  I  1)*C2(  12)  *C3(  13)  *F(  11,12,13) 
IF  (13  -1)117,  116,,  117 

116  TSM=0.0 

117  TSM  =  TSM  +  T  +  TS 

103  Z( JU,JV,JW) =TSM 

107  (MARK  =  MARK  +1 
JJW  =  JW  -  1 

PRINT  5,  JJW,  U,V,W,  Z(JU,  JV,JW) 

WW(JW)  =  JJW 

IF  (XMODF  (MARK,  35))  109,108,109     ^ 

108  IF  (NO  134,  15,  133 

134  PRINT  8 
PRINT  6 
GO  TO  135 

133  PRINT  8 
PRINT  2 

135  PRINT  4 

GO  TO  165  p 

109  IF  (XMODF  (MARK  ,  5))  165,110,165 

110  PRINT  7  „ 
1-65  CONTINUE 

104  CONTINUE 

IF  (TRY)   72,44,44 
44   GO  TO  3 
72  'IF  (ND)   163,163,73 

7  3  CALL  DRAW( 101 ,Z ,WW ,0 ,0 , LABEL , ITITLE,0,10.,0,0,0»0»6»10,1,LAST) 
GO  TO   163 

END 
END 
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